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ВВЕДЕНИЕ 


В данных методических указаниях содержится описание 
четырнадцати работ, помогающих исследовать методы решения задач 
численного анализа, с тридцатью различными вариантами заданий. 

При выполнении работы необходимо предварительно 
ознакомиться с соответствующим методом и его алгоритмом, отчетливо 
представлять себе исходные данные задачи, процесс решения и 
ожидаемые результаты. Далее должна быть составлена блок-схема, по 
которой пишется программа на языке высокого уровня. В качестве 
языка программирования может использоваться любой язык, например 
Паскаль. Отладка программы на ЭВМ с анализом полученных 
результатов составляют заключительную стадию работы. Возможна 
проверка полученного результата при помощи специализированных 
математических программ, использующих язык прикладной задачи. 

При составлении отчета в нем, после титульного листа, 
указывается вариант задания, используемые вычислительные формулы, 
блок-схема алгоритма, текст программы на языке высокого уровня, 


результаты и выводы по используемым методам поиска решения. 
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ЛАБОРАТОРНАЯ РАБОТА № 1 


Решение нелинейных уравнений 


Нелинейным уравнением с одним неизвестным называется 
уравнение вида: 


Ро) =0 (1.1) 


где Р(х) - непрерывная,  дифференцируемая функция, 
определенная на некотором бесконечном или конечном интервале 
а<х<Ь. 

Корнем уравнения (1.1) называется такое значение х, равное & 
которое обращает уравнение в тождество: 


ғ =0 (1.2) 


Найти корни уравнения (1.1) точно удается лишь в частных 
случаях. Кроме того, часто уравнения содержат коэффициенты, 
известные лишь приблизительно, и следовательно, сама задача о точном 
определении корней уравнения теряет смысл. Поэтому разработаны 
методы численного решения уравнения вида (1.1), которые позволяют 
отыскать приближенные значения корней этого уравнения. 

При численном решении уравнения (нахождение корней) в общем 
случае необходимо проделать три процедуры: 

а) установить количество корней и область значений аргумента, 

которой принадлежат эти корни; 

б) выделить участки изолированности корней, то есть такие 

участки, внутри которых расположен только один корень; 

в) сузить изолированный участок до заданных размеров, то есть 

найти корень с заданной точностью. 


При пользовании большинством методов решения уравнений 
предполагается, что первые две процедуры выполнены, то есть известен 
участок изолированности корня. Графически корень нелинейного 
уравнения изображается точкой пересечения графика функции }(х и 
оси ОХ (рис.1). 

Условием существования корня уравнения (1.1) на промежутке 
(а, Б| значений аргумента х для непрерывной функции | (х) является 
неравенство: 


Ла Ў) 50 (1.3) 


Геометрически оно обозначает, что функция (х) на концах 
промежутка [а, В] имеет значения разных знаков. 
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Если (х) на промежутке [а, В] не только непрерывна, но еще и 
монотонна (то есть ее производная не изменяет знака), то тогда, при 
выполнении условия (1.3), на отрезке |а, Р| содержится только один 
корень уравнения (1.1). 





Рис. 1. Корень нелинейного уравнения. 


1.2. Метод проб 


Метод проб является одним из простейших численных методов 
поиска корней. Суть его состоит в том, что некоторый заданный отрезок 
значений аргумента [о1, рі| разбивается на небольшое число 
подотрезков, и для каждого из них проверяется выполнение условия 
(1.3). Тот из подотрезков, для которого условие (1.3) выполняется, 
содержит корень, поэтому его разбивают на более мелкие подотрезки. 
Процедура повторяется до тех пор, пока не будет найден подотрезок 
достаточно малой длины (соизмеримый с заданной погрешностью), 
содержащий корень. 

В качестве расширения метода можно предложить схему, в 
которой после нахождения корня аналогично производится проверка 
оставшейся части интервала [а, В]. Такая модификация метода хороша 
тем, что позволяет найти все корни на заданном промежутке. 
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1.3. Метод половинного деления 


Начало 


Описание Г(х) 
определение 
а, 6, є 





Рис. 2. Блок-схема алгоритма метода половинного деления. 
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Для применения этого метода необходимо, чтобы функция на 
отрезке [а, Б| содержала только один корень. При нахождении этого 
корня отрезок [а, В] делится пополам, то есть выбирается начальное 
приближение х= (а+0)/2, затем выбирается тот из отрезков |а, х| или 
[х, р], на концах которого функция } (х) имеет противоположенные знаки. 
Новый отрезок снова делится пополам и так далее, до тех пор, пока 
длина отрезка, на концах которого функция имеет противоположенные 
знаки, не будет меньше заданного числа =. Блок-схема метода приведена 
на рис. 2. 


1.4. Метод Ньютона 


Метод Ньютона (метод касательных) состоит в том, что 
приближенным значением корня Х считается точка пересечения 


касательной к кривой } (х) и оси ОХ (рис. 3). 





Рис. 3. Корень уравнения, определенный методом Ньютона. 


Для применения метода необходимо: 

а) чтобы на отрезке [а, В] производная ў ' (х) была всюду отлична от 
нуля; 

б) чтобы функция {(х) была на [а, Б| монотонна (производная 


сохраняла зн ак) 5 
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в) в качестве начального приближения следует выбирать такую 
точку хо, для которой выполняется условие / (хо) - Ё" (хо) > 0 (то 


есть точку, в которой кривизна и функция имеют одинаковые 
знаки). 


Рабочая формула метода 


ЛО) 
Хи = Ха = (1.9) 
Р(х) 
позволяет найти последовательность Хо, Хі, ж, ..., хь, сходящуюся к 


точному значению корня 6 уравнения (1.1). Критерием окончания 
вычислительного процесса является выполнение условия: 


[хамы -Хп| <=. 


Алгоритм метода 


1. Выбрать начальное приближение хо, принадлежащее 
промежутку [а, Ы|, положить х = хо. 
ХО) 


2. Вычислить А = 





709 
3. Если |Л| < є, то корень С = х-А, вычисления окончены. 
4. Положить х = х-А и вернуться к шагу 2. 


1.5. Метод итераций 
Метод итераций (последовательных приближений) предполагает, 
что уравнение (1.1) преобразовано к виду: 
х= ф (х). (1.5) 


Последовательные приближения корня хо, хі, ж, ... задаются 
рабочей формулой: 


Хп = ф (хп) прип = 0, 1, 2,... (1.6) 


Последовательность Хо, Хі, Х2, ..., Хп, ... сходится к точному 


значению корня & если для любого х є |а, Р| выполнено условие 
сходимости метода: 


ІФ (| «1 (1.7) 
Алгоритм метода 


1. Выбрать начальное приближение х є |а, БІ. 
2. Вычислить Л = ф (х). 
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3. Если |х-А| < є, где є - заданная точность, то & = Л и вычисления 
окончены. 


4. Положить х = Ли перейти к шагу 2. 





Рис. 4. Последовательное нахождение корня методом итераций. 


Вопрос о преобразовании уравнения (1.1) к виду (1.5) решается в 
каждом отдельном случае по своему. В общем случае можно 
рекомендовать некоторый универсальный, но довольно громоздкий 
прием. Если ищется корень уравнения (х) = 0 на отрезке |а, Б, то 
необходимо найти число 


= зр ГО) 





; (1.8) 


хе| а, 


равное значению точной верхней границьт модуля производной функции 
709 на отрезке [а, Б]. Тогда ф (х) можно выбрать в виде: 


ф(х) =х-2709); а 


при этом условие сходимости (1 7) выполняется автоматически. 


Варианты заданий приведены в табл. 1. Для всех заданий принять 
точность вычислений = = 10-5. 
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Таблица 1. Варианты заданий к лабораторной работе № 1. 
№ Отрезок, Метод № Отрезок, Метод 
Тех Уравнение содержащи | численного ат Уравнение содержащи | численного 

Р. й корень решения Р. й корень решения 
половинного а: половинного 

[0;0,85] итерация х+с08(х052+2) =0 [0,5;1] итераций 

8 +5103, З+еіп3,6х > 
1055 1: 25 1 1 половинного о 2 й половинного 
половинного половинного 
= а &] 2 -1: 
сов 25 Е а _ аа а Е ый 
1 
[0;1] Ньютона 26 хівх-- зб [0,2;1] Ньютона 
деления 2 деления 
2 

ох . половинного д. половинного 
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ЛАБОРАТОРНАЯ РАБОТА № 2 


Отделение и уточнение корней нелинейных ин Описание / (х). Задание 
отрезка, количества участков 


уравнении определения разбиения, точности, 
определение шага 


Данная работа аналогична работе №1 с той разницей, что здесь 
Определение 





задан ориентировочный участок функции, где могу быть один или границ 

подотрезка 
несколько корней. Поэтому необходимо сначала отделить и определить 
достаточно малый участок, где есть корень, методом проб, а затем 
уточнить значение корня быстросходящимся методом. Так как корней на Корень есть? 
заданном участке может быть несколько, например, как показано на Да 
рис. 5, то после уточнения корня необходимо проверять наличие корней а 

ьвод 

и далее до конца участка. контрольных / 77777777 Вывод концов отрезка, на 


значений котором находится корень 


При выполнении задания для метода проб отрезок предлагается 


делить на 20 частей. Из-за сложности анализа функций на сходимость 
Уточнение Метод половинного 


для уточнения корнеи использовать метод половинного деления с корня ПИ деления или какой-либо 


заданной точностью. другой 


оны И для контроля значения 
функции в найденной точке 


Изменение 
номера 
подотрезка 





Рис. 5. Отделение корней нелинейного уравнения. 





Все под- 
отрезки? 


Блок-схема алгоритма приведена на рис. 6, варианты заданий в 


таблице 2. 


Рис. 6. Блок-схема алгоритма отделения и уточнения корней 
уравнения. 
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Таблица 2. Варианты заданий к лабораторной работе № 2. 


№ Заданная 
Уравнение Отрезок 
вар. точность є 
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2 
ЗОВ ОН СИ 
роза реж 
соѕх+ Ме ря Ух =0 [2; 8] 
Е р | 
нано |а оно 
З аа 
мем 
ВЕ 
п о 
ШЕ 38 
И 





10922 - Л+ю (1+) =0 
вро 
сем ЕВ МЕССИ СИИ 
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ЛАБОРАТОРНАЯ РАБОТА № 3 


Вычисление определенных интегралов с 
заданным шагом интегрирования 


Численные методы нахождения значений определенного интеграла 
р 
І- Ї Го) ах используются в тех случаях, когда неизвестна 
а 


первообразная для функции (х) или ее вычисление сопряжено с 
трудностями. Большинство методов основано на замене 
подынтегральной функции аппроксимирующей функцией более простого 
вида и последующим интегрированием отрезков этой более простой 


функции. 
1. Метод прямоугольников. 


Этот метод основан на аппроксимации кусочно-постоянными 
функциями. Рабочая формула метода: 


п-1 
1=һУ ЛО кА) 

> (3.1) 
а. а 

п 





Здесь п-число отрезков разбиения. 

В случає, если: 
Д=0 - получаем рабочую формулу левых прямоугольников; 
Л=й/2 - формулу центральных прямоугольников; 
Д=й - формулу правых прямоугольников. 


2. Метод трапеций. 


В этом случае аппроксимация производится кусочно-линейной 
функцией в соответствии с рис. 7. Рабочая формула метода: 


п-1 
І «Годі Гуд +25700) 
і=1 


а а х= а+іћ 
п 


(3.2) 





Где п – число отрезков разбиения. 
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Рис. 7. Замена значения интеграла суммой площадей трапеций. 


3. Метод Симпсона. 


В этом случае подынтегральная функция аппроксимируется 
отрезками парабол. Рабочая формула метода: 


п-2 п-1 
ЕЕ У Л) +4 У, С 





і-20) із 10) (3.3) 
р-а 3 
ћ = ‚ хізачій 
п 
Здесь п - число отрезков разбиения (обязательно четное). 
Формулу можно привести к виду 
ћ п-1 ; 
1 = а ло ОЕ» 3-р |700} (3.4) 
1=1 


Блок-схема алгоритма примера приведена на рис. 8, варианты 
заданий в табл. 3. Вычислить приведенный в варианте интеграл 
методами центральных прямоугольников, трапеций и Симпсона при 


одинаковом заданном количестве отрезков разбиения. 
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Начало 


Описание Кх) 

Определение 

значений а, Б, 
п, й 





Вычисление 
суммы $ 
(циклический 
процесс) 


Интеграл = 
п/З [Ка) +60) + $] 







Вывод 
значений 
интеграла 


Рис. 8. Блок-схема вычисления интеграла методом Симпсона. 


Таблица 3. Варианты задания к лабораторной работе № 3. 


Количество частей 


разбиения 
Шх 
=== 30 
1 ХУ1+ 0 х 
л ІЗ 
Ї (5? х+ се? х) ах 54 
л16 
36 





Вычисление определенных интегралов с заданным шагом интегрирования 19 


Количество частей 


сеет разбиения 


4р2 
Уя 


1 х 


122 
| Мех —1 ах 
0 


1 

[ хе" ѕіп х ах 

0 
2 
[хѕћх ах 48 
0 


зы 1 


(ав 
0м9+ х2 


1 РА 


х 


= 
| х агсів х ах 
0 


У х 


3 
[га + х) ах 
1 


1 1 


А 
01+ 3х+ 2х2 


ах 





20 


Га 


хагсіє х 


0+ х2 


ах 





Лабораторная работа № З 


Количество частей 
разбиения 


160 


240 


22 


48 


22 


96 


60 


52 


176 


36 


Вычисление определенных интегралов с заданным шагом интегрирования 


1,5 
Ї віп х Ш (є х) ах 


3/4 


р сн 
0 (х+1)Ух2 +1 


1 





1 
| 5 ах 
о (3Ззшх+2с0$х) 


2 


1х4 +3х2 +2 


21 


Количество частей 
разбиения 
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ЛАБОРАТОРНАЯ РАБОТА № 4 


Вычисление определенных интегралов с 
заданной точностью 


Обычно при вычислении интегралов заранее известна точность, с 
которой необходимо произвести расчет. Количество же участков 
разбиения зависит от свойств подынтегральной функции и выбирается 
намного больше, чем это необходимо для обеспечения заданной 
точности. 

Чтобы не делать лишних вычислений (а с ростом вычислений 
увеличивается и аппаратная ошибка, связанная с конечной длиной 
числа, представленного в ЭВМ), используют алгоритм с удвоением 
отрезков разбиения и сравнением интегралов, вычисленных с п 


отрезками разбиения и 2 п. 


Алгоритм метода 


1. Вычисление интеграла с двумя отрезками разбиения. 

2. Увеличение отрезков разбиения в 2 раза. 

3. Вычисление интеграла с новым количеством отрезков. 

4. Сравнение предыдущего и вычисленного значения интеграла. 
Если разность значений больше заданной точности, то перейти 


на шаг 2, иначе конец алгоритма. 


Если вычисление интеграла (шаг 3) производить по алгоритму, 
приведенному в работе №3, то вычисление одних и тех же значений 
функции будет проводиться многократно, что не рационально. 
Желательно при вычислении интеграла с новым количеством отрезков 
использовать результаты предыдущих вычислений. Это достигается 


следующим образом . 
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На рис. 9 представлена схема вычисления интеграла по 
приведенному выше алгоритму. Очевидно, что значения (хо), и Ё (хп) 
постоянны для интеграла с любым количеством отрезков разбиения, 
поэтому их можно вычислить один раз в начале программы. Значения 
промежуточных результатов (х! для двух отрезков разбиения) будут 
использоваться как четные слагаемые для следующего значения 
интеграла (х22 для четырех отрезков разбиения). Таким образом, остается 


вычислить только сумму для нечетных значений функции. 


Блок-схема алгоритма для данного примера с использованием 
метода Симпсона приведена на рис. 10. Формулы интегралов для 
вариантов заданий приведены в табл. З к лабораторной работе №3, 


точность вычислений для всех вариантов принять 10-5. 


Хх (х5) 





Рис. 9. Вычисление интеграла с двукратным увеличением 
количества отрезков разбиения. 
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Начало 


адание 
подынтегральной 
функции, границ 
интегрирования, 
количества точек 
разбиения (2 


Вычисление 
2=Кх)+Кх,) и шага 
р, четное значение 
суммы $ч=0, нечет- 
ное $„=/((хо+х,„)/2) 


Вычисление 
старого значения 
интеграла 
Ч/З (2+4$н) 


Увеличение 
количества точек 
разбиения в два 

раза и вычисление 
шага 
интегрирования 


Определение 


четной суммы 
54=54+5 н 
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Вычисление 
нечетной суммы 
(в цикле) 


п-1 
У Р(х) 


1=[2) 


Вычисление нового 
интеграла по 
формуле Симпсона 






Заданная 
точность 
достигнута? 






Нет 


Новое значение 


интеграла 
становится старым 





Конец 


Рис. 10. Блок-схема алгоритма вычисления интеграла с заданной 


точностью. 
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ЛАБОРАТОРНАЯ РАБОТА № 5 


Поиск корней сложной функции с 
использованием подпрограмм 


Часто при решении задач численного анализа используются 
кусочно-непрерывные функции, функции, описанные рядом, интегралы 
и так далее, которые нельзя описать в одном операторе. Поэтому 
необходимо использовать подпрограммы, в частности подпрограммы- 
функции. Примером данной задачи является поиск корней уравнения, 
заданного интегралом. 


Пример. 


Найти корень уравнения (а) =0, соответствующего функции 
согласно варианту задания, методом половинного деления на указанном 
отрезке с заданной точностью. 

Здесь вычисление (а) необходимо оформить в виде 
подпрограммы-функции с аргументами а и є (точность вычисления 
интеграла). Интеграл вычисляется методом Симпсона с заданной 
точностью. Начальная точность вычисления интеграла 0,2, затем 
увеличивается в 2 раза на каждом шаге по мере уменьшения отрезка, на 
котором находится корень. Когда корень найден, значение функции 
выводится с точностью 105-5. 

Программа должна состоять из основной программы и 
подпрограммы-функции. 

Основная программа составляется по блок-схеме рис. 2, 
приведенной в лабораторной работе №1, а описание функции 
оформляется в виде подпрограммы. Так же необходимо в подпрограмме 
предусмотреть увеличение точности вычисления интеграла (не путать с 
точностью вычисления корня) при каждом последующем приближении. 
Значения аи (а) (контрольные точки) выводить в каждом цикле, чтобы 
определить правильность работы программы. 

Подпрограмма-функция составляется по блок-схеме рис. 10, 
приведенной в лабораторной работе 4. Входными параметрами будут 
являться верхняя граница интегрирования и заданная точность. Все 
остальные значения, в том числе и подынтегральная функция, могут 
задаваться локально внутри подпрограммы. Варианты заданий 
приведены в таблице 4. 
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Таблица 4. Варианты заданий к лабораторной работе №5. 


Точность 


Отрезок, 
нахожде 


содержащий 
корень 


а ш0,5х+е05х 
га) = | з 


01 50х+1,5 





ах 


а 
Га) = | х0 (п 0,5х+ х) ах 
01 


Раз [е —0,2х?) ах 


_ Я 0,45х+е04* 
бл зіх-15 


ах 





а 
Га) = | хіп 0,6х + х) ах 
0,1 


уо) = [ея — 0,42) ах 


0,1 





а 
(а) = [х0 (а 0,7х+ х) ах 
0,1 


Газ | (е7 —0,6х?) ах 
0,1 


Е | 10,35 х + е°.2* 


01 зшх+1,5 


ах 
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Отрезок, 
содержащий 
корень 


а 
Јо) = | хо п 0,8х+ х) ах 
01 


Раз [кеа — 0,8х?) ах 


а ш0,Зх+ех 
Га) = | 2802 


01 $Шх+15 


ах 





а 
Јо) = [х0 0,9х + х) ах 
01 


п 0,25х + е0! 
01  8іпх 15 





а 
Јо) = | Зп 0х + х) ах 
0,1 





а 
Јо) = | х0 (а х+ х ах 
0,1 
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Точность 
нахожде 





28 


пауз | є" -1Ах)) ах 
0,1 





а т 0,15х не Зх 
ах 


ее 


01  8іпх 15 


а 
Јо) = [хо (в 120+ х) ах 
0,1 


(о) = (ее —1,6х?) ах 


а 01+ е0 


Га) = | 


01 $Шх+15 





ах 


а 
Га) = [х0 (а 13+ х) ах 
0,1 


(ед = [кее —1,8х2) ах 


а 110,05х+е0°* 
Ра) = | а 


01 50х+1,5 





ах 


а 
Га) = | х0 (п 14х+ х) ах 
0,1 


Газ | (е2 —2.0х?) ах 
0,1 
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Отрезок, 
содержащий 
корень 


Точность 
нахожде 





Решение систем линейных алгебраических уравнений 29 


ЛАБОРАТОРНАЯ РАБОТА № 6 


Решение систем линейных алгебраических 
уравнений 


Системой линейных алгебраических уравнений п-го порядка 
называется система 


а11м + ато +... + ах, = р 
ауд + 422 В.п = (6 1) 


ат + аро +... Фр = 0, 


Решением этой системы является значение неизвестных хі, ж, ..., 
Хһ. Коэффициенты матрицы а; и правые части р; предполагаются 
заданными. Условием существования и единственности системы (6.1) 
является неравенство нулю определителя матрицы 


деё а#0 (6.2) 


В дальнейшем предполагается, что это условие всегда выполнено. 


1. Метод Гаусса (схема единственного деления) 


Идея метода Гаусса основана на последовательном исключении 
неизвестных и преобразовании системы (6.1) к специфическому виду 
(6.3) системы с треугольной матрицей (прямой ход метода): 

м + ах + аж +...+ ах Е 50 


хо + аб) 3 +... + ах, = 2 


зара КО 25. (6.3) 
(п. р(п=р) 

п-1,пХп Е 1 
в Е 


Хи— +1 


После этого осуществляется последовательное определение 
неизвестных хи, Хп, ..., Хі (обратный ход метода). 

Метод Гаусса относится к группе так называемых «точных» 
методов, теоретически он позволяет с помощью конечного числа 
действий вычислить точный результат. Однако при практических 
расчетах возникают вычислительные погрешности, которые искажают 
результат тем сильнее, чем выше порядок системы. Условием сходимости 
метода является отличие от нуля определителя матрицы системы. При 
программировании метода и отладке программы целесообразно по 
окончании прямого хода вывести для контроля треугольную матрицу. 
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Алгоритм метода 


1. Положить $ = 1. 


а . 
$ 5 
2. Для ]= 5 +1, $+2, ..., п положить ау о, 
ас 
Ь 55 
3. Положить В; =—®. 
485 
. Положить і = 5+1 
‚Для ј= 5+1, 5+2 ,..., п положить а; = а 45 ау; 


. Положить р; = Бі - Бе аз; Із і+ 1. 


. Если 1< п, перейти к шагу 5. 


очо пт 


. Положить $=$+1; если $ < п -1, то перейти к шагу 2, в противном 
случае к шахту 9. 


9. Вычислить х, =—" 





аһ 


10. Положить $ = 5-1 
п 
11. Вычислить х, -ф,- Уашх 


К=5+1 
12. Если 852 2, то перейти к шагу 10. 
2. Метод простой итерации 


Идея метода состоит в том, чтобы, задавшись некоторыми 


0) 


значениями неизвестных м (начальным приближением) : С 


помощью вычислительной процедуры получить «уточненные» значения 


неизвестных ж зх Р 


Это уточнение продолжается до тех пор, пока 
значения хм... хи не будут вычислены с заданной точностью. Для 


применения метода система (6.1) должна быть преобразована к виду: 





м = В! з дохо — 13% нин 
хо = Во - 091 - 223Х3 —--.— Фәх 





(6.4) 





хи = В, ал 7 22 — --- 7 0 п-1Хп-1 
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Рабочая формула метода: 





к 2 ї 
х Тов У аух() 
ј=1 
ар. : 
іж 
аз = 4а; Ј (6.5) 
0, із) 
Ь. 
В; = — 
дії 


для 1=1,2,...,п;]=1,2,..., п. 


Условиє сходимости метода од «1, где |< |- норма матрицы аў. 


Критерий окончания вычислительного процесса 


«ен _ Ю 
1 


: «є, 


тах |х 
і 


где = - заданная точность вычислений. 

Методы простой итерации относятся к группе «приближенных» 
методов. Это значит, что они дают за конечное число действий лишь 
приближенный результат, однако этот результат может удовлетворять 
пользователя своей точностью. Критерием сходимости методов так же 


является условие е1 <1. Под нормой матрицы можно понимать любое из 


следующих чисел: 


|А = па а 6.6) 


(6.7) 


в 
| 


(6.8) 


е8 
| 
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При программировании метода полезно организовать выдачу на 
экран всех последовательных приближений.  Задавая значения 
коэффициентов із необходимо убедиться в выполнении условия 


сходимости. В качестве начального приближения, если нет иных 








соображений, можно выбрать нулевые значения хі = 0 = хз=... = Ха =0 


или значения правых частей хі = В+. 
Алгоритм метода 


1. Зададимся значениями хі, ж, ..., хп. 


2. Для #1, 2, ..., п вычислим значения у; = 8; — р Е і: 


З. Найти величину 6 - наибольшую из разностей | х:- у;| для і- 1, 
2л ПО 

4. Если б < є „то у; - искомые значения, и вычисления можно 
окончить. 

5. Для 1- 1,2,..., п присвоить х: значение у.. 


6. Перейти к шагу 2. 


Решить систему уравнений, с коэффициентами, заданными в 
таблице 5, методами Гаусса и простой итерации. Во втором случае 
предварительно определить норму матрицы. Если не будет выполняться 
условие сходимости, то путем линейных комбинаций в системе 


уравнений привести ее к сходящемуся виду. 
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Таблица 5. Варианты заданий к лабораторной работе №6. 


Матрица коэффициентов 
Свободные члены 
чо 
— 


_ 
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ЛАБОРАТОРНАЯ РАБОТА № 7 


Определение собственных векторов и 
собственных значений матрицы 


Собственными значениями матрицы называются корни Л, А», .. 
Ап характеристического уравнения п-го порядка 


её |А-ЛЕ|, (7.1) 


р) 


где А - матрица пхп. 
Е - единичная матрица порядка п. 

Собственные значения матрицы могут быть упорядоченными по 
величине ЛМ > 22» ...> А, ставится задача отыскания наибольшего 
собственного значения Дах. 

Собственным вектором матрицы А, соответствующим і - тому 


собственному значению А: называется решение х; системы 


В силу того, что система однородна, собственный в ектор 


х; определяется с точностью до произвольного постоянного множителя. 


Поскольку любой произвольный п - мерный вектор а может быть 


разложен по собственным векторам х; 
= п = 
аз Усіх, (7.3) 
і=1 
то при подстановке вектора а в систему получим 


п п п 
Е =: = К да 
Аа= АЎ сх, з Усх, = | сх, 6, Р 
і=1 із! і=2 1 


й 


= пар п А, а 5 = п А, => 
Аа= А сах + Ус, х, = 4 ах +Ус, Б 
ї=2 Я, 1=2 1 
з 3 
сли з Еа А, Кес 
А а= А сх + Ус — | х; и так далее. 
1=2 
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В общем случае 


ро ДО у 

А а= А1 Сх] + |2) хі. (7.4) 
22 “А 

жк 


1 


Если |4121401> 14| >... > 4, ,то величины -»0 при К> о то 








есть при достаточно больших К 
Аќа = Д ем 


или Аа = ле хуя А“ ад. 


Таким образом, выбрав произвольный вектор а и, умножая его 
последовательно на матрицу А К раз, получим, что А! является 
коэффициентом пропорциональности в соотношении 


Аа = ДАКа, (7.5) 
причем точность этого равенства повышается с ростом К. 
м ЦЕ - 
Тогда 4 = = ›где | - норма вектора х. 
|А а] 


Алгоритм метода 


1. Выберем произвольный вектор ио = (1,0,...,0). 


2. Вычислим вектор и! = Аи) и первое приближение собственного 





значения да = ‚ где |а | = 





и 


3. Вычислим вектор из = Аи и второе приближение собственного 

















и 
значения Ж -ТПОМ 
и 
2 
4. Если |? д? |<= то о — приближенное значение первого 
собственного значения, из - приближенное значение 


собственного вектора. Если |А - 2 | >= то положить и зи», 
АЭ = 49 и перейти к шагу 3. 


Варианты заданий приведены в таблице 5 к лабораторной 
работе 6. Заданная точность вычислений є = 103. 
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ЛАБОРАТОРНАЯ РАБОТА № 8 


Интерполяция сеточных функций полиномами 


Простейшая задача интерполирования заключается в следующем. 
На отрезке [а, В] задано п+1 точек хо, хи, ..., Хи, которые называются 


узлами интерполяции, и значения некоторой функции } (х) в этих точках: 


Р(х) = уо, Г (ха) = чл, ..., (ха) = уһ, (8.1) 


Требуется построить функцию  Е|х) (интерполирующая функция), 
принадлежашую известному классу и принимающую в узлах 
интерполяции те же значения, что и /Д/(х. При использовании 
интерполяционной формулы Лагранжа функция Е (х) является 


полиномом Ѓл (х) степени не выше п, который выражается следующей 





формулой: 
п (х= хох ху о хи ЗО У Ху) 
І = Я ДЕО а БА ом 2 8.2 
9 б 2 (х; ді хох; = хр)...(х; = хх; = х41)...( Хі Е. хи) | | 
или 
п п ХХ, 
Г.о) = Уу, П---- (8.3) 
1=0 ј=00— Ху 
)кі 


Таким образом, для решения данной задачи необходимо 
организовать тройной циклический процесс, два цикла которого 


предназначены для вычисления одной точки полинома [Гм (х). 


Используя интерполяционную формулу Лагранжа, вывести в виде 
таблицы значения функции ј (х) для точек х с заданным шагом, если 
известны значения функции / (х) в узлах хі. 

Блок-схема алгоритма для данного примера приведена на рис. 11, 


варианты заданий в табл. 6. 


Интерполяция сеточных функций полиномами 


Начало 


Определение 
массивов, 
начальные 
занесения 








Начальные 
занесения 


Начальные 
занесения 


Вычисление 
полинома 





Вывод 
текущей 
точки 


Рис. 11. Блок-схема алгоритма 


Перебор точек, в которых 
необходимо определить 
значение функции 


Вычисление суммы 


Вычисление произведения 
с проверкой 


жі 


интерполирования функции. 
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Таблица 6. ик заданий к лабораторной работе № 8. 
Отрезок 
— О А і ана 
у Е Е О ВЕСИ ри ре БАН 
хі Я 1,5 | 2,4 | 3,8 6,8 | 7,1 | 8,2 
а аа |=. 
хі 1,5 | 2,1 | 2,9 | 3,1 | 3,4 | 3,8 | 4,6 | 5,2 
3 хі 1,7 6 7,4 | 10,5 [12,8| 15 | 16,1| 19,8 4 18 1 
у | 5,7 |6,4 | 6,5 | 4,2 | 2,1 5 9,8 116,6 
хі |0,06 | 0,15 | 0,221 0,28] 0,48 |0,57 | 0,7 | 0,89 
хі 1,1 | 1,3 [1,6 [2,1 | 2,2 | 2,6 | 2,9 З 
х |8,11 |8,44 | 8,7 |9,01 |9,55 |9,66 | 9,9 |10,2 10 | 01 
ш | 4,2 | 7,1 | 5,4 4 6,7 | 7,3 | 8,6 | 9,2 й 
хі |0,1|0,3|0,75| 1 11,1111,37 |1,72| 1,9 
хі |0,99(1,13(|1,2511,39| 1,5 | 1,62 | 1,78 | 1,88 
а Р 3 я о | м |і 
х |2,3 | 2,8 | 4,1 | 4,5 | 5,2 | 5,9 7 7,9 
ОНЕ 
хі |0,77 | 0,9 | 0,94] 0,9610,99| 1 1,17| 1,3 
ВЕЕР 
Хі О |0,25|0,35| 0,5 | 0,75 | 0,95 | 1,2 | 1,5 
Хі -1 0 0,7 | 1,2 | 1,9 4 4,3 | 4,9 
[563 а аз 2 [88 9 [лао 
хі | 0,5 | 1,2 | 17 | 2,2 | 3,1 | 3,5 | 3,9 | 4,1 
Хі 1,3 | 1,9 | 2,8 | 3,7 | 4,2 | 5,1 | 5,8 | 6,9 
[393888 | я [881882 | $ [62 
м | 0,3 | 0,9 | 1,1 | 1,8 | 2,1 | 2,8 | 3,4 | 3,9 
о ооо | 3 [88 з 38 о | а | а 
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№ Отрезок 
Ў Значения 1 интерполирован 


вар иа 

х |0,1 | 1,5 | 1,8 | 2,2 | 2,9 | 3,4 | 4,5 5 

хі 1,2 | 1,8 | 2,6 | 3,1 | 3,6 | 4,5 | 5,6 | 6,2 
ЕЕ: 


Пе › |5 
оаа 
ЗЕ ЕЗ ЕЕ ЕЕ ЕН ЕЕ С 
ра ре реа р а ро р фе 
ЗП ЕЕ ЕВ ЕЕ ЕС З КАС 
ЕЕ ЕЕ В ЕЕ ЕН ЕВ С СОС 
о 0 ЕЕ ЕЕ КЕ СО СС 
ааа оја 
оаВааававааоо 
КЕНЕ ЕЕ ЕЕ ЕЕ ЕС 
СЕ ЕЕ ЕН ЕЛЕ ЕВ ЕВ З 
новее 
п аа аз] З ва 
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ЛАБОРАТОРНАЯ РАБОТА № 9 


Среднеквадратичное приближение сеточных 
функций 


п 
Дана сеточная функция у (х), заданная на системе узлов {х;},_| 


п 
своими значениями в узлах Ут: Известно, что по системе п узлов 


можно построить интерполяционный полином порядка п-1 , который 
будет принимать во всех узлах в точности предписанные значения. 
Однако, во многих практически важных случаях необходимо построить 
полином степени меньше, чем п-1, приближающий данную сеточную 
функцию «в среднем». В этом случае можно сформулировать задачу 
отыскания такого полинома порядка т < п-1: 


т 
Р. (х)= У а х! =а,+ах+а,х? +азх+..+а.х" | 9.) 
1=0 


при котором «невязка», то есть суммарное отличие значений полинома в 
узлах от значений сеточной функции, было бы наименыпим. Если 
записать выражение для «невязки» в виде 


п 


2 
(ааа: [р (х,)- у, | , (9.2) 
1=1 
то условиями для определения козффициентов полинома ао будут 


необходимые условия минимума функции многих переменных 
5 (ао, ал, ..., ап). 


Таким образом, коэффициенты {а ии могут быть найдены из 


решения системы уравнений 


08. =0, 7=0,1,..,т (9.3) 
да, 


Для получения рабочих формул метода продифференцируем 


выражение для «невязки» 5 = У [р, (х)- у | по ај: 


і=1 
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25 2 











Ш а 2 
Б У Ре) - 3 = >, 10) - УИ" = 
ба; да; г 1 да; 
М й | . (9.4) 
п А 
= Ув, ори 
р дај 
Поскольку Р, = У; арх! то ВЕ) =х/ и условие =0 
ж 7 б 
приводит к уравнению: 
п і . 
2УТРи(ху) - 315 =0, 7)-0,1,..т. (9.5) 
1=1 


Преобразуем выражение суммы 
п і п т і і т п ко) п 7 
У, [2„0)- Ух = Х| Ў ам - У |х = У а Ў - ЗУ (9.6) 
ыі ЕО зар чаї а 
Следовательно, неизвестные коэффициенты ао, а, ..., ат МОЖНО 


найти, решив систему из т + 1 линейного уравнения 


п В п Е п А п і 

+1 +т 

[ай а За ун Ха )- У умі і (9.7) 
1=1 і=1 і=1 і=1 

Алгоритм метода 

1. Зададим узловые значения функции у; и узлы х для і= 1 + п. 


2. Вычислим матрицу й размером (т +1)х(т +1): 


п А а 
ар = Ух |, дая і ј= 1+ (т +1). 
К=1 


3. Вычислим правые части системы 0:4 =: 


п 
1—1 А : 
= Ууру , для = 1+ (т +1). 
К=1 
4. Решив систему уравнений методом Гаусса, получим ее решение: 


а, ј= 1+ (т+1). 
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5. Вычислим значения найденного полинома в узловых точках 
т-1 


ЕЕ 1-1 а. ра 
Б. = Хар; , і-ізп 
7-1 


и вычислим «невязки» в узлах б; = В _У. 


і 


Варианты заданий приведены в табл. 7. Значения найденного 


полинома степени т искать только в узловых точках. 


Таблица 7. Варианты заданий к лабораторной работе № 9. 
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Значения 1 


8,44 9,01 | 9,55 | 9,66 | 9 
7,1 4 6 8 
0,43 1 
6,7 5 
1,13 
1,7 
2,8 
9 
0,9 


Значения і 


40 43 
5 -3 
1,12 | 1,25 
16 -29 


28 32 
12 
1 
37 
8 11 13 
38 72 61 
2,9 
72 


8,11 7 й 5 9 
55 4,2 5,4 7 | 73 ‚6 
0,96 0,21 0,75 1,11 | 1,37 | 1,72 
47 6,9 1,3 |7 9,4 | 7,6 2 

А 1,25 З 1,62 | 1,78 

4 1,4 9 1 -3,3 | -9,6 
_ 3 4,1 2 |59 | 7 

2 7,2 -2,9 8 


1 


6 

2 
1, 
2 

-1,2 | -0,9 | -0,5 

-0,6 | -2,2 


2 
5,85 


4 
83 


8 

З 

1,39 | 1,5 

2 3 | -1,6 
4,5 2 
4,3 | -0,1 


2 


7,3 


23 
99 
0,9 
20 
5 
46 
0 
16 


» фо + ! 


Е 
ое | + Ф Ф 
оо | б м 


ро 
оо 


4 
5 


1 
ро 
№) 


2 


на 
е1 


2 
6 
2 
6 
7 
6 


6 
2 
4 
4 
9 
9 
4 
5 
6 


с 





4,7 

10,5 | 12,8 

4,2 2,1 
0,15 | 0,22 | 0,28 | 0,48 

-2,8 | 2,2 0,5 


0 


1,2 
1 
70 
2,8 
15 
1,3 
0,3 
1,5 
2,4 
1,7 
5,7 

06 
0,6 
1,1 

1 


м 
00 


15 
2 
1, 


№ 
с 


9 
28 
4,1 
6,4 
1,5 
2,2 
2,1 
5 
6 
6,4 
2,7 
1,3 
8 





1 2,1 2,2 
-4,4 | 3,6 | -8,7 
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ЛАБОРАТОРНАЯ РАБОТА № 10 


Шаговые методы интегрирования 
дифференциальных уравнений первого 
порядка 


Дано дифференциальное уравнение первого порядка у = ѓ (х, у (х)) 
с начальным условием у (хо) = уо. Необходимо вычислить сеточное 
представление у: частного решения у= у(х), удовлетворяющее 
дифференциальному уравнению и начальному условию. 

Численные методы решения задачи Коши делятся на шаговые и 
разностные. К шаговым относятся методы Эйлера, Рунге-Кута, Гилла. К 
разностным - Адамса, Милна, Штермера, Гаусса и другие. 


1. Метод Эйлера 





Рис. 12. Решение дифференциального уравнения методом Эйлера. 


Решение задачи в конечном итоге сводится к интегрированию 
обыкновенного дифференциального уравнения 


а 
ууу) (10.1) 
ах 


с начальным условием 


у (хо) = уо. (10.2) 
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Если принять / (х, у (х)) = сопѕі (метод левых прямоугольников), то 
Убх) = Уб) + В Рх, УС) (10.3) 
— рабочая формула метода Эйлера. 
Алгоритм метода 
1. Задаем ћ, начальные значения хо, Џо, 1. 
2. Вычисляем уін. 
3. Изменяем іи х. 


4. Если найдены все точки, то конец, иначе перейти к шагу 2. 


2. Метод Эйлера-Коши 


При интегрировании дифференциального уравнения (10.1) по 


формуле трапеций, получаем уравнение 





Уж +) = уб) +5 но) НУВ, УЛ (0104) 


Выражение ц(х; + П) из правой части находится по методу Эйлера 
(10.3). Таким образом, рабочие формулы метода Эйлера-Коши имеют 
вид: 

ж 


Ура = УА Р(х, ур) 


я 10.5 
Уна = У «За У ТИМ уг) зам 








Алгоритм метода Эйлера-Коши аналогичен алгоритму метода 
Эйлера. 


3. Модифицированный метод Эйлера 


В данном методе при интегрировании уравнения (10.1) 


применятся формула центральных прямоугольников: 


убогу (у(х 5) (10.6) 


48 Лабораторная работа № 10 


Неизвестная з +5] находится также по методу Эйлера (10.3). 
Рабочие формулы модифицированного метода Эйлера: 
А ћ 
Уну = У: +5 (ху) 
(10.7) 
й * 
Ум = УС Д х, Фо Уну 


Алгоритм аналогичен алгоритму метода Эйлера. 
4. Метод Рунге-Кутта четвертого порядка 


Из широкого класса методов типа Рунге-Кутта наиболее широко 
используется метод Рунге-Кутта четвертого порядка. 
При разложении у (х + П) в ряд Тейлора в окрестности точки хі до 


четвертого порядка малости получаем: 





Е ей: 
убх +) = Ус) ВУС) + у) + уто) +217 (цу) (10.8) 
2 б 24 
Решение ищется в виде Ц(х: З И) = у(х) + Ді, 
где Л: = а К: + а2Ко: + азКзі З адКаї. 
Итоговая рабочая формула: 
К = ВРО, У;) 
ћ К; 
Юю; "С 90 ш) 
ћ К; 
Кз; "е З +2) 
Ка = В У(х; +, у; + К) 
Ка; +25; + 203; + Кд; (10.9) 
Уц = У; + І 21 Зі 4і 





6 


Алгоритм аналогичен методу Эйлера. Варианты заданий 


приведены в таблице 8. 
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Таблица 8. Варианты заданий к лабораторной работе № 10. 


№ Дифференциальное |Начальные Отрезок Метод 
вар интегри интегриров 


уравнение условия 
рования ания 


ия 
со8х 
2}, . Модифиц. 
(+х2)у’+1=0 и0)=1 | [6:05] | 0,025 | "Зара 
' 9-х Е Й Эйлера- 
у +2у = 2е  соѕх и(0)=1 [0;0,5] | 0,025 
' д Е . Модифиц. 
6 у= у =5іпх+со52х 0,1 
' = 25 У Модифиц. 
а й і й аа 2 Г й а 
Е убу =2е“* /(0)=1,433 | [0;0,5] | 0,025 
Е РУ 7 г Модифиц. 
подо Е 4 Эйлера- 
4 =: = _ Я Модифиц. 
" ДЕН - зай . Эйлера- 
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Отрезок Метод 
интегри интегриров 
рования 


Зу+2у=х 2 Зх [0;1] Рунге-Кутта 
У 


ва 


Дифференциальное |Начальные 
уравнение условия 


(0)= 
(0)= 
2у' – = у= 8х [1;1,5] 
| оне СОС 
Г и вы фе 
ом о Да До 
а нер Це |з 
р око нара |а 
х 
у +4у = Зѕіпх+ѕіп 2х РОС 
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ЛАБОРАТОРНАЯ РАБОТА № 11 


Интегрирование дифференциального 
уравнения разностным методом Адамса 


Введем сетку значений аргумента хі = хо + іП, где ћ = соп5і - шаг 
сетки, 1= 0, 1, 2, ... Сеточное представление частного решения у: = у (хі). 
Для его определения существуют шаговые методы, в которых очередное 


сеточное значение определяется по данному: 


узі = У (уй), (11.1) 
и разностные методы, в которых очередное сеточное значение 


определяется по нескольким предыдущим 


уза = И (у, ум, уго, ...). (11.2) 

Формула разностного метода Адамса может быть получена, если 
проинтегрировать исходное дифференциальное уравнение на отрезке 
[х хні| и подынтегральную функцию в правой части уравнения 
заменить интерполяционным полиномом Ньютона для интерполирования 


в конце таблицы. 





х1 х1 АЕ 42 і 
ун оз Јоу) урн коду 25 

х; ее ћ р 

і б (11.3) 
х(х—)(х— х1) + Ес со 2) +. ах 


Ограничиваясь  полиномом третьего порядка и проведя 
интегрирование полинома, получим рабочую формулу, в которой 
значение конечных разностей заменены их выражениями через 


значения в узлах сетки: 


ун = 21557: - Аа] (11.4) 
Р = Е 
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При использовании разностного метода Адамса проблемой 
является получение первых, так называемых «стартовых» точек, в нашем 
примере это уо, и1, у2, уз. Они могут быть получены с помощью шаговых 
методов точности не ниже третьего порядка, например методом Рунге- 


Кутта четвертого порядка. 
Алгоритм метода. 


1. Определяем стартовые точки у, іл, 1р, уз. Присваиваем 








значения А = уо, В=ш, С = 10, "” = уз. 
. Вычисляем ја = }(ж, А), јв= Р (хи, В), с = (кю, С, Љ= р(х, Р). 


3. Вычисляем новые значения сеточной функции 


№ 


ћ 
УР - 59 /с +37 в -9ХА| 


4. Если хз + А> хкон, то конец вычислений, в противном случае 


перейдем к шагу 5. 








5. Переприсваиваем хз = хз + А, ја = в, в= јо, јс= р, = Л (ж, у), 
"Р = у. 
6. Переход к шагу 3. 


Дифференциальные уравнения к вариантам заданий приведены в 


таблице 8 лабораторной работы №10. 
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ЛАБОРАТОРНАЯ РАБОТА № 12 


Поиск экстремумов многомерной функции 
методом градиента 


Методы оптимизации (поиска экстремумов) очень широко 
используются на практике. Существует большое разнообразие численных 
методов оптимизации, различающихся по размерности решаемой задачи 
(одно- и многомерные), по наличию ограничений (безусловные и 
условные), по способу формирования шага (градиентные, 
безградиентные, случайного поиска) и по другим признакам. Как 
правило, численные методы отыскания экстремума состоят в построении 


УЮ 
последовательности векторов Х 5 при поиске минимума 


удовлетворяющих условию: 


р (12.1) 


Обычно ищется минимум, отсюда и название - методы спуска, а к 
максимуму перейти достаточно просто. 

В данных методических указаниях рассматриваются только 
градиентные методы без ограничений. В этих методах элементы 


(К) 
последовательности їх вычисляются по формуле: 


(+0 з2(Ю 509 
Х = Х +арР й (12.2) 


К =0, 1, 2, ... п 


500 
где Р - направление спуска, ак- длина шага в этом направлении. 
Направлением спуска является антиградиент: 


роб - х 9) (12.3) 


а длина шага в различных методах определяется по-разному. 

Метод градиента в чистом виде формирует шаг по переменным 
как функцию от градиента заданной функции {в текущей точке поиска. 
Простейший алгоритм поиска минимума для двумерной функции 
записывается в скалярном виде в соответствии с формулой: 


4 
а = 5 

а О З (12.4) 
Упа == й—— 
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Величина рабочего шага в направлении антиградиента зависит 
как от самого градиента, так и от коэффициента пропорциональности Й, 
с помощью которого можно управлять эффективностью метода. 
Существуют разновидности метода как с постоянным, так и с 
переменным шагом. Последние здесь не рассматриваются. 

Для оценки частных производных в общем случае используются 


разные методы с двумя вариантами: 


1. Алгоритм с центральной пробой 


аї - а Ен хо х , Хінно х,) | 
ах. 5; 


1 





2. Алгоритм с парными пробами 


аї Д Р(х, а РЕАЛ х) РР х,) 


2 





ах; 2.8; 


1 


где д: - пробный шаг по Кой переменной, выбираемый достаточно малым 
для разностной оценки производной. 
Условием окончания поиска минимума может являться условие, 


когда величина модуля градиента менее заданной ТОЧНОСТИ 5. 
Пример 


Требуется найти минимум функции 
Лоьу)= (х2 + ›?) 


с заданными начальными координатами х= 0.5; у=0.5 и 
точностью вычислений є =10-3. 

Зададимся значениями пробного шага д=10-3 и коэффициента 
шага А 0.3. Обычно данный коэффициент приходится подбирать 
вручную. Так, при ћ =0.1 минимум находится за 32 шага, а при ћ =1 
процесс не сходится. Поэтому надо вводить вторую проверку окончания 
цикла: по количеству приближений. 

Если выводить промежуточные результаты вычислений, то 


получим следующую таблицу: 
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Таблица 9. Пример вывода результатов к работе № 12. 



































№ х у ар/ ах а ау [атаа Досу; | Лх,у) 

1 0.5000 0.5000 0.8776 0.8776 1.2411 0.4794 
2 0.2367 0.2367 0.4705 0.4705 0.6654 0.1118 
З 0.0956 0.0956 0.1911 0.1911 0.2703 0.0183 
4 0.0382 0.0382 0.0765 0.0765 0.1082 0.0029 
5 0.0153 0.0153 0.0306 0.0306 0.0433 0.0005 
6 0.0061 0.0061 0.0122 0.0122 0.0173 0.0001 
7 0.0024 0.0024 0.0049 0.0049 0.0069 0.0000 
8 0.0010 0.0010 0.0020 0.0020 0.0028 0.0000 
9 0.0004 0.0004 0.0008 0.0008 0.0011 0.0000 
10 0.0002 0.0002 0.0003 0.0003 0.0004 0.0000 





























Использовать для четных вариантов точность &=10-3, алгоритм 
оценки частных производных - парные пробы. Для нечетных вариантов 
точность є =10-2, алгоритм оценки частных производных - центральная 


проба. Варианты заданий приведены в таблице 9. 
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Таблица 10. Варианты заданий к лабораторной работе № 12. 
Начальное 
№ Функция приближение 
хо Чо 
1 (оа ноя) уніа) |ноазадногв у? (1.25 | -1.95 
0.75 е 
чіп [0352+0.2.(у-0.32}} + 
2 -2.5 5.2 
+50з(0.47-(+0.25]+0.34-(у—0.21} +0232 
зав (0.22-(х-0.1}+02.(у—0.32} 
З -2 5 
+ 0.43 х2+1.97 (0.34) –0.45 ху 
4 | 20.085х2+0.12+0.022%у + ві (0.22х+0.8у+1.3)+0.43х2 | -3 | З 
5 2 0.033х2+0.028у2-51п(0.38х-1.12)_ 024 у2+0.69ху 4 11 
6 е9 х2+у+іп(1.8х-1.12у) е02 у2+х+со1.44х+1.85у) -2.5 0.2 
2 2 
0.43 х +0.32 у" - 0.45 хуч 
7 б і З зі 
+98 [032+ 605(0.34-(у-0.21] 




















бек кетанов иномаркой ФУАТ ТОМІ ИКАНА Лабораторная работа №8 









































№ Р Начальное 
иззи приближение 
0.42 х2+0.37 у2+ 
8 2 2 -3 0 
зи 005 (х-0Л8 ов Озон і; | 
2 2% 
9 „№023х +0.26 у +іп(0.8х-—1.1)+-соѕ(1.1у+0.75х) 0.7] ху -4 1 
10 (на (ноу 1.2--023х24021 у2+0.75ху 21 | 0.05 
0.82 -е 
зіп (0.3.(5-0.13}+0.2-(у+0.22} + 
11 -2 -1.2 
+3.1-с0(0.88 22+0.78-(у—0.12)? )+0.39 х2+0.46 у? 
зп(0.67-(:—0.85)2+1.2-(-13) + 
12 0 5 
+2.73-(х+0.25)2 + 3.07 (30.77) 1.79 ху 
13|] 20.1222+0.1у?-0.0820у 4 віп (| 2х 0.35у -1.3)- 035 -2.5 | -0.5 
14| г 0.093х2-іп(0.89х-—1.02)+соз(1.55у+0.2х)  ().88 у? +0.45ху 0.8 -5 
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№ ення Начальное 
2 У приближение 
15 д 0.82х2-у+віп(2.8х+2.2у) +. 077 у2-х+со$1.4х+1.8у) 0.5 2.6 
2; 2 
0.13 х" +0.22 у" - 0.45 хуч 
16 -2 1 
+58 [0.6-(х=0.39} )+-соѕ 0.8-(у-0.2)) 

1.72 х2+1.07 у2+ 
17 2 2 2 0 

зн [048 (54084085 зовуть р | 

2 26, 
18 „№9575 +0.63 у нап (1.25х+1.3)+со1.1у+0.9х) 0.25 ху -4 0 
19 А || (ному? уніа)? |537+6+08752+009 2 12 | 19 
е 

хап ((—0.73)20.33 з?) 
20 1.5 1.3 

+с0з[0.51-(х-0.79}+1.35-(у-1.38} }+0.78 х2-+0.91 у? 

виа [0.5-(5+0.25}}—1.13-(5+0.88 + 
21 0 З 
+2.02 .(х+0.89)2 + 2.77 -(3+0.51)+1.66ху 

22| 20033:7+0.02857-00320у 4 5п(1.09х+0.75у-1)-04у2 | 8 | 10 
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Начальное 
№ Куса приближение 
23 о 0-132—51п(0.8х—1.1)+со(2.15+0.75х) 12 у2+0.71ху 4 су 
24 б 0.23х2+у+ѕіп(2.8х+3.2у+1) е 018 у2+х+со3.4х+3.7у) 0.5 0.5 

0.15 х2+0.2 у2+0.28 хуч 
25 -5 -5 
+511(0.06-(х—0.18)? + соѕ[0.098-(у-031)2) 
1.1х2+1.23 у2+ 
26 9, 2 -1 -1 
зві ооо ов чол О | 
27 4(0.57х2+0.63 у? +оп(1.67х+0.3)+соѕ(1.985+0.66х) 0.3 ху 2 -4 
Р т 
28 | сем у-02)? |2.97+3+0.155:2+0097 у2 09 | оз 
е 
зав [3.45-(х-1.03}2+0.1.(у-0.67} + 
29 1.5 1.3 
нсо*((х-0.09}+2.3.(у—1.08} +2.27 х2+2.72 у? 
зав (1-7. х—0.15}2-1.4-(у—0.5 ) | 
И (х-0.15)?-1.4-(у-0.5) И 
+ 2.8. (х+0.3}2 +3.7-(у+0.7Р+2.7ху 
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ЛАБОРАТОРНАЯ РАБОТА № 13 


Поиск экстремумов многомерной функции 
методом тяжелого шарика 


Данный метод базируется на аналогии с движением материального 
шарика по наклонной поверхности. Шарик будет стремиться занять 
нижнее положение, то есть точку минимума, со скоростью в зависимости 
от крутизны наклона, его массы и вязкости среды. Эти параметры 
влияют на характер движения шарика. 

Для двумерной функции траектория поиска в дискретном 


варианте описывается в соответствии с (13.1). 











аў 

Хр х; а(х; хц) ћ у 
і рен ДА, (13.1) 

=у-а(и а) 

Урни Уау У ду, 


При а=0 метод превращается в обыкновенный градиентный, при 
а=1 процесс поиска не затухает. То есть эффективность метода зависит 
не только от ћ и начальных условий, но и от 0< а <1. Первый шаг 
делается обычным методом градиента, так как значение предыдущей 
точки еще неизвестно. 

Варианты заданий приведены в таблице 10 к лабораторной работе 
№12. Использовать для нечетных вариантов точность є =10-3, алгоритм 
оценки частных производных - парные пробы. Для четных вариантов 
точность є =10-2, алгоритм оценки частных производных - центральная 
проба. Так же необходимо ограничить количество шагов на уровне 30- 


100. 
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ЛАБОРАТОРНАЯ РАБОТА № 14 


Поиск экстремумов многомерной функции 
методом наискорейшего спуска 


Основным недостатком градиентного метода является 
необходимость частого вычисления производных от заданной функции. 


В методе наискорейшего спуска длина шага ак определяется из условия: 
(к (6 (к зл 
тех нок за) (14.1) 
а20 


то есть на каждом шаге находится оптимальное значение а, мини- 
мизирующее одномерную функцию по направлению антиградиента: 
один раз вычисляется градиент и ищется минимум по этому 
направлению. 

Для двумерной функции Ё(х, у) ак определяют следующим 
образом. После определения градиента минимизируемой функции 
делается шаг в направлении, обратном градиенту и пропорциональный 


ему (чем меньше крутизна функции, тем меньше и шаг): 


* А 
то зх, – 90} 


З Е (14.2) 
Уп+1 = Уп — 90у 





Если значение целевой функции при этом уменьшилось по 
сравнению с исходным 
* * 
ГО Уһ) > ХО Уһ) › 
то делается очередной рабочий шаг в том же направлении: 


ж А р 
аа захо ідо Їх 


Е | | За (14.3) 
Упчі = уп — 0 Лу 
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Если же после очередного шага значение целевой функции 
увеличилось по сравнению с предыдущим значением, то движение 


прекращается: 


ак = @-1, 


# 
Хи Хата 


Уп+1 = Уп За 





и заново определяется значение градиента. Перед началом следующего 
цикла проверяется ТОЧНОСТЬ нахождения минимума, которую 


приближенно можно оценить по величине градиента: 


Сн =? (14.4) 


Определенные трудности представляет выбор ао, от которого 
зависит сходимость метода. Если ао мало, то процесс будет сходиться 
медленно, в противном случае может оказаться, что при первом же 
рабочем шаге в очередном приближении значение функции увеличится, 
и новое приближение будет равно старому. Поэтому необходимо 
предусмотреть проверку: ак должно быть больше нуля. 

Кроме этого, количество приближений п должно быть меньше 


заранее заданного числа М. 


Блок-схема алгоритма для данного метода приведена на рис. 13. 








При выполнении работы считать = =0,1; №М=20; ао =0,1. Частные 
производные определить аналитически. Если процесс расходится, 
необходимо считать с двойной точностью. 


Варианты заданий представлены в табл. 11. 


Начало 


Определение 
Қх,у), (ху) 


Їу(%,у), Хо, Уо 
П, Ч, &, 


Определение 
градиента в 
точке х и у 

п п 












Достигнута 
точность 
вычислений? 


Нет 


Определение 
старой 


величины 
функции 





Очередной 
рабочий шаг, 
определение 

новой величины 


функции 





Функция 
бывает? 


Нет 


Новая 
функция 
становится 
старой 
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ппре- 
высило 
предел? 









Да 


Нет 


Определение 
ак 











Определение 
НОВОЙ ТОЧКИ Изменение 
приближения ао 


Хинт И У пої: 





Вывод а к , ко- 
ординат точки, 
функции и 
градиента 


Рис. 13. Блок-схема алгоритма метода наискорейшего спуска. 
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Таблица 11. Варианты заданий к лабораторной работе № 14. 


Минимизировать функцию 


Лоьу)=2х-1Зу+ 


Минимизировать функцию 


Убу) =16х+е199:7+026)° 


2 2 
РО, у) =17х+ Олуже256 +0,27у 








2 2 
ЈО, у) =3х-1,2 Оу) =18х+0,2у +289: +028 
х,у =3х- 1, уче 











2 2 
РО, у) =19х+0,3у+е324х +0,23у 











ГО, у) =4х-Шу+е 


ЛОу)-5х- ручне) 250 +01532 


Р(х, у) = 20х+0 дез? +0,3у2 








уу) = 21х+0,5у+е2 "02197 
РО, у) = 6х 0,9у+ ©0536? +01 6у2 





ЧА Р(х,у) = 22 х е502х? +0,32у2 
ХУ =/х- > уче 


Е 2 
2 

| ЈО, у) =8х-—0,7у + е0,64х +0, 

Е да 











2 2 
РО, у) = 24х 2529х+0,34у 





Ле, у) =9х-0,6у+е 











гоху) 255 20,35у2 





2 2 
Р(х,у) = 10 х - 0,5у+е0 94 +0,2у 


ГО, у) = 6х+ужеб25* +0,36у? 


2 2 
РО, у)= Пх-0,4у+ех +0217 


2 2 Р(х,у) = 27 х У Ыу цебто +0,37у2 
РО, у) =12х-0,3у+ е!2 1х7 +0,22у 





2 2 
Оу) = 28 х+Ь2у 26.98” +0,38у 





2 2 
ГО, у) = 13 х-0,2у Ср га +0,23у 











ГО, у) = 29х 27.2957 н0,39у2 


Оо, у) = З0х+1.4у-не А +0,4у? 





2 2 
(х,у) = 14 х –0/у +ег69х +0,24у 


2 2 
ГО, у) =15 х+е!9% +0,25у 


= о, у) = 23 х 24.84%? +0,33у2 


2 
З 
4 
6 
й 
10 
11 
12 
13 
14 
5 





